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We study dynamics of rolling tachyon and Abelian gauge field on unstable D- 
branes, of which effective action is given by Born-Infeld type nonlocal action. Pos- 
sible cosmological evolutions are also discussed. In the Einstein frame of string 
cosmology, every expanding flat universe is proven to be decelerating. 



1 Introduction 

Recent analysis of string theory l ^ l ^ l suggests that classical dynamics of 
tachyon on an unstable D-brane is described by a scalar Born-Infeld type 
action with a runaway potential El. Characteristic of such tachyon effective 
theory is absence of perturbative tachyonic modes around the potential mini- 
mum and existence of spatially homogeneous rolling tachyon configuration El. 
In the background of the rolling tachyon, Abelian g auge f ield on the D-brane 
cannot support perturbative spectrum at late-time * ' ' V In this review, we 
summarize firstly bosonic effective action both in the bulk and on the D-brane, 
together with shapes of the tachyon potential. In section 3, we introduce 
Carrolian limit, and then discuss nonexistence of perturbative spectra of the 
tachyon and the Abelian gauge field in the potential minimum, of which dy- 
namics are governed by nonlocal D-brane action. For the tachyon, there exists 
the classical rolling tachyon solution connecting smoothly unstable maximum 
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of the potential and stable minimum at infinity. 

Though the situation seems not to be matured yet, string cosmology must 
be an intriguing subject to be tackle d at every step of progress in both string 
theory and observational cosmology Recently, inspired by string theory 
D-branes and heterotic M theory, brane cosmology attracted much attention. 
Brane cosmology assumes that our universe starts out with branes embedded 
in the higher dimensional spacetime, either stable or unstable. I n such context, 
cosmology of the rolling tachyon must be an attractive subject I ' ' ' ' '. 
Topics included inflation, dark matter, cosmological perturbation, and reheat- 
ing despite stringent difficulty in the simplest versions of this theory. A nice 
aspect for the cosmology is nothing but the absence of open string excitations 
so that classical analysis of late-time based on homogeneity and isotropy may 
lead to a solid prediction. In section 4, we briefly describe system of gravity 
and the tachyon. The obtained universe in the simplest set predicts that ex- 
panding uni vers e is accelerating at onset stage but turns to decelerating phase 
at late-time Finally, the scale factor approaches a constant. Our main 
concern is string cosmology involving graviton, dilaton, and the tachyon. In 
the string frame, we obtain flat space solutions at both the early epochs and 
late-times in closed form. When transformed to the Einstein frame, analysis 
shows that every expanding universe of the graviton-tachyon-dilaton system 
should be decelerating irrespective of the specific shape of the tachyon poten- 
tial. Specific form of the scale factor is proportional to (time) 1 / 2 . 

We conclude with a summary of this paper. Sometimes, same notations 
are used in different sections for convenience. 



2 Setup 

Effective action of bosonic sector of D3- brane system of tension T3 is described 
by the following action in string frame ISID 

S = i / dix ^~9 e~ 2 * (R + 49^$ - -^H^H^r - 2a) 

-T 3 J tfx e -*V(T)yJ- det(g„ v + d^TdJT + + F M „) , (1) 

where cosmological constant A is assumed to be vanished in what follows. 
The (quadratic) action in the first line of Eq. JQ is bulk action of graviton g^, 
dilaton $, and antisymmetric tensor field B^ u of rank 2 of which field strength 
is defined by H^ p = d^B up + d v B pil + d p B^. The action induced on 
the D3-brane is Born-Infeld type nonlocal action given in the second line of 



"Other effective actions different from Eq. can also be considered. For example boundary 
string field theory (BSFT) action for a non-BPS D-brane was also proposed and studiedly 



2 



Eq. |Q, which involves kinetic terms of tachyon T and Abelian gauge field 
satisfying F^ u — d^A v — dyA^. Note that the tachyon T and the Abelian 

gauge field A^ were rescaled as T — > \/2na'T and — > 2na' ' A^ to absorb 

27ra' in front of the kinetic terms. 

Since we will deal with two classical top ics in this paper, i.e., one is 

string cosmology including the rolling tachyon H^l and the other is evolution 

of electromagnetic waves in the background of the rolling tachyon ^ let us 

firstly read equations of motion from the action Q 

R„ v - \g^R = -A<? M „ - 2 [V M V„$ + g^ (V Q $V Q $ - V Q V a $)] 



+k T 3 e V(T) I — 9n„V9 x ), ( 2 ) 

4(-V Al V^$ + V Al $V"$) = R-2A-j-H a ^H a ^-K 2 T 3 e* ^(TV^T 1 *, 

(3) 

Vf,H^ x = 2V (ti M 71!/A , (4) 

C?KC a PV u.TV v X a a - XWTV V C?K 
Cffi^T ^ ^ - C^V^TV^ 

i^f^-C^TV,T)=0, (5) 



V dT \ ^ 

where X^ = + + iv„ + 8,/rduT, X = det(X^), W is the cofactor 
of X^ v , and is the symmetric part of C". 

The tachyon potential V(T) multiplied in front of the Born-Infeld type 
action Q is known to measure varying tension of the unstable D-branes. 
Though there seem no consensus and no exact computation on the form of 
tachyon potential V(T), it should have its maximum, max(V(T)) = 1, and 
minimum, min(V(T)) = 0. Some known computations from string theories 
are as follows^! For bosonic string theory, we have 

V{T) = / ? 6XP i 1 ~ f ° r SmaU T but T ~ 1 (6) 
[ exp (— -j) for large T 

and, for superstring theory, we have 

for small T but T > 
--4= I tor large L 
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V(T) 




String theory calculation © _ tells us that favorable form of the tachyon 
potential supports V ~ e~ T for small T and V ~ e _T for large T. Since 
the potential V(T) must be a smooth function in the inte rme diate region, 
connecting two asymptotic expressions, a typical example is 1 "^ 

^ = ( 1+ I) KP (-|)- ,8) 

In this paper, considering convenience in analytic computation, we adopt 

cosh (^) 

which is also consistent with the results of string theory ©-0 (See Fig.QJ. 



3 Tachyon and Gauge Field 

Among various bosonic degrees of freedom, dynamics of the tachyon T and 
the gauge field is governed by nonlocal kinetic terms in the brane action 
(|T|). Before discussing their dynamics directly, let us consider a massless 
nonintcracting real scalar field cj> described by Lorentz-invariant linear wave 
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equation 



( 




(10) 



If we take two singular limits known as the In6n6-Wigncr limits, infinite ve- 
locity limit, v — ► oo, goes to Galileo's limit of action at a distance described 
by an elliptic equation : 



which we encounter frequently in Newtonian dynamics. On the other hand, 
the opposite velocity dominated limit, v — > 0, called as Carrolian limit is 
described by an ordinary differential equation : 



A natural implication from the above discussion is possible absence of 
plane waves in the Carrolian limit of corresponding field theory. Instead time 
evolution of homogeneous configuration may be dominated. Though both 
limits are singular in the context of group theory, we will use the speed v as 
a continuous parameter of range (0 < v < 1) as we have done in the Galileo's 
limit where its value is taken from from one to infinity. 

In this section, we will briefly review rolling tachyon in flat spacetime 
with disappearance of perturbative tachyonic spectrum at the bottom of the 
potential V(T) and discuss singular evolution of perturbative electromagnetic 
field in background of the rolling tachyon. 

3. 1 Rolling tachyon in flat spacetime 

Let us begin this subsection with discussing physics of unstable D-branes 
in the simplest setup. If we consider solely motion of the tachyon in flat 
spacetime, it is depicted by El 



where specific form of the tachyon potential V(T) is chosen by Eq. ©. Note 
that g _1 det(g /J j / + d^Td^T) = 1 + d^Td^T up to the quadratic terms and it 
holds in an exact form as far as homogeneous configurations are concerned 
such as the cosmology described in terms of Robertson- Walker metric of our 
interest. Equation of motion is 



V 2 </> = 



(11) 




(13) 




(14) 
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where fj = -(1 + d^Td^T) and 

d^Td v T 

v = v + W ^ = ?r ______ (15) 

Energy-momentum tensor 7),„ is given by 

Tfj, u = -TsViT)^^^^. (16) 

Since T = is unique stationary point of the tachyon potential V(T), we 
perturb the action i|13|l up to quadratic terms to read perturbative spectrum. 
For small T, we have 



•S ^ / * T,, I </'./■ f -ld»Td»T + ±sT 2 j . (171 



where the positive cosmological constant at T — is subtracted in the action. 
As expected, the sign in front of mass term is negative so it is tachyonic 
on top of the potential V(T = 0). For large T, the tachyon potential is 
approximated as V(T) « 2e _T / T °. Introducing a new small order parameter 
such as <f> = 2v / 2Tbe~ T / 2T ° , we obtain a quadratic action with positive mass 
term 



T 3 / d 



j4„ 



(18) 



Though this seems to imply existence of a perturbative tachyonic spectrum 
of mass l/\/2To, it is an artifact of quadratic approximation since T — > oo is 
not a stationary point of the tachyon potential as shown in the Fig. ^ To be 
specific, let us consider a plane wave solution <f) ~ e ik i* x and then we have 
d^T ~ —2iTok fl . Inserting these into the equation of motion 114|l . we get the 
following dispersion relation 

1 dl » V 0. (19) 



y/1 - AT Q 2 k 2 dT 

Note that d(\n V) / dT\T^oc = —To ^ 0, the dispersion relation does not have 
any nontrivial solution but can only be satisfied for infinite k 2 . The mass-shell 
condition we derived through a quadratic approximation <|18|) was reproduced 
by a Taylor expansion of Eq. I|19|) up to quadratic terms for small k 2 . 

Absence of perturbative tachyon spectrum does not lead to nonexistence 
of classical configurations. We may take into account an intriguing lesson from 
the Carrolian limit (|12|) of the wave equation. Suppose that we have a ho- 
mogeneous tachyon configuration with d{F = and consider its classical time 
evolution T(t). Then only diagonal components of the energy- momentum ten- 
sor (llfit survive and the equation of motion (|14|l is integrated to conservation 



G 



of energy px = Tqq, i.e., Hamiltonian density H.t is 



H t =Pt=T 3 V( - T \ = T 3 e, (20) 
\/\-T 2 



where e = V(Ti)/Jl - Tf > 1 and T 4 , T; are the initial values of T and T. 
It is analogous to the aforementioned Carollian limit Ijl2|l . Constancy of the 
Hamiltonian density (121 III forces T — * 1 for any finite e as T — > oo. For the 
tachyon potential ©, we have a solution of rolling tachyon 

T(t) = T Q sinh- 1 Le*/^ - , (21) 



where a = ~ 



sinh(T t /To) + ^smh 2 (Ti/T ) + b 



and 6 = (e 2 - l)/e 2 . At late- 



time regime, the obtained tachyon solution (|21|1 is approximated as expected, 
i.e., T(t) w t irrespective of the initial values of T and T. Note that the 
pressure pr = T\/3 is negative and approaches zero as time elapses, which 
coincides with vanishing Lagrangian limit 



p T = C T = -T 3 V{T) Vl-T 2 '-^ 0. (22) 
From equation of state px = wtPt, we read 

w T = -(l-T 2 ) = ~ V{T f )2 <0, (23) 
where air = — 1/e 2 at T = and wj- — > at T — + oo. 

3. ,3 Gauge field and rolling tachyon 

When we turn off antisymmetric tensor field of second rank B^ v in the Born- 
Infeld type action in Eq. JTJ, effective bosonic action of unstable D3-brane 
system is given by 

S = -T 3 J d 4 x V(T)y/- det^ + d^Td v T + F^) . (24) 

To proceed, we introduce a few notations. We first define 

= Vtu , + d^Td u T + F„„, (25) 
X = det(X M „), (26) 

and barred covariant field strength tensor 
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with Eq. (|15fl . Then we have contravariant barred field strength tensor F^ v 
and its dual field strength F* 

= fj^fj^Fcp, F; v = e J^lp a(3 = e -t^lfj a ^f s F j5 , (28) 

where e pi , Q/3 = ^rj e^p with e i 2 3 = 1. 

In terms of barred quantities Eq. I|26[) is computed as 



i + \h v F^ - 1 (f; v f^y 



2 ' 16 

Then equations of motion for the tachyon T and the gauge field are 



d, 



V 



CP* d v T 



V 



-X 



dT 



eg* ] = 0. 



(29) 

(30) 
(31) 



Here Cg" and are symmetric and asymmetric part, respectively, of the 
cof actor, 

C» v = fj (fT + F^ + ff a ff 1 ff v Fl i pF* b + f\» a ^KpF; 5 F Sv ) , (32) 
namely, 

Cf = fj(fT + fj^f^K^s), (33) 
= fj{F^ + r a n Pl F* p F; s F Su ). (34) 

Energy-momentum tensor is given by 

T 3 V(T) 



-r]^X + C w {d v TdPT + F U P)] 



(35) 



where C M „ = r]^ a ri u pC a ^ . The diagonal components can be identified as 
density and pressure, 



P 



T 3 V(T) 
T 3 V(T) 



X + Co^fd^T + F^) 



(36) 



-3X + C^{fd»T + F ") + C^{d^Td v T + F^) . (37) 



3V-X 

From equation of state p = wp, we read 

1 AX -C^{8^Td u T + F^) 



X + C 0fi (Td»T + F^) 



(38) 
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Thus the trace of the total energy-momentum tensor is given by 



T 3 V(T) 

V^x 



l-i + W^TO'T) -!- j (1 + ^d p Td p T)6» a 



d^Td a T Fu V F' 



(39) 

The equations of motion (|30|l and l|31|) admit homogeneous solutions T{t) 
and A^{t) which satisfy 



do 
do 



V 



-X 



-.T 



0. 



V 



-X 



E l = 0, 



(40) 



where X = — 1 + T 2 + E 2 . Solving these equations we find that V/ \J — X = 
constant and E = constant. The first equation is nothing but the statement 
of conservation of energy. 

When E = 0, the solution reduces to that in the theory with pure tachyon 
field given in the subsection 3.1. When E ^ 0, one can simply read the solution 
from Eq. (2D by rescaling e 2 -> e 2 (l - E 2 ) and T -> TqI^/Y^W- in the 
exponents of Eq. (f?T|l . 

Now we consider the propagation of the gauge field A p in the rolling 
tachyon background (|21|l . For this purpose we linearize the equation of motion 
(|31|l in this background. Written in terms of electric and magnetic fields, 
Ampere's law is modified as 



V x B 







V 



vVi - f 2 dt \Vi-f 2 



E 



(41) 



Combined with the Bianchi identity, the wave wave equations of the electric 
field E and the magnetic field B become 



1 



d 



V dB 



vVi-t 2 dt V Vi-f 2 9t 



d_ 

at 



V\/l - T 2 dt V V 7 ! - f 2 



-V 2 B = 0, 



V 2 E = 0. 



(42) 



(43) 



Note that these equations are not identical. Since the Hamiltonian density 
Ht (|20jl is a constant of motion, the wave equations and reduce to 



8 2 B 



e 2 9 2 E 
V 2 ~W 



V 2 B = 0, 



^l^l^-V^E^O. 



d 
dt 



V 2 dt 2 
2 \ 9E 



e 

V 2 



(44) 
(45) 
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One can also compute the energy density l|36|) and pressure l|37[) as are given 
by 

T 3 V(T) ( E 2 



Pc 



Tn 



3 (T« 



T, 



00 F„ 



7ii|.F„»,=o) 



2 v/T^T 2 
T 3 V(T) 



1 -T 2 



6 v 7 ! -T 2 



(1 - T 2 )B< 



Similarly, we read w from Eq. I|38|) 



Pcm 
Pcm 



Fern ^ j>2\ 



The trace of the energy-momentum tensor l|39[) becomes 



rp )X 

em fi 



rpfl 



t 3 y(T) 



rp2 



E 



(46) 
(47) 

(48) 

(49) 



The speed of wave propagation is altered by y(T(i)) and thus time- 
varying. One might expect from Eq. (|44|l that the propagation speed roughly 
follows V(t), but this is not necessarily so as we will see below. Consider- 
ing the shape of the potential in Eqs. © and (J7J), we see that it is the light 
speed on the top of the potential and decreases monotonically to zero as the 
rolling of the tachyon proceeds. For illustrative purpose, we take the tachyon 
potential JHJ and the solution l|21|l . Inserting the solution into the tachyon 
potential, we have an expression of V as a function of time t 

V(T(t))= 1 (50) 

y/l+ (ae*/ T ° - s e ~* /T °) 

As shown in Fig.0 V(T(t)) changes in the scale of a few Tq's from the light 
speed to zero, which is nothing but elapsed time from the onset of the tachyon 
rolling to the late-time. Since the electric field and the magnetic field obey 
different wave equations, they show different behaviors in propagation. 

To investigate in detail the propagation of electromagnetic wave in the 
rolling tachyon background, let us try to find a plane-wave solution of the 
wave equations with a time- varying propagation speed, (|44|l and (|45|) . We 
consider the plane-wave propagating along the x-axis, k = fcx, and make an 
ansatz B = B z (t,x)z with 

B z (t,x) = B{ty k[x ~ WB{t) \ (51) 

where B(t) and Wb(£) are real. We assume that V^(0) « 1 and Wb(0) ~ 1 so 
that the ansatz l|51|) gives the normal plane-wave initially. Since E and B are 
related through the Ampere's law l|41|) and Faraday's law 



VxE+ — 
ot 



0. 



(52) 
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5 10 15 20 

t/T 



Figure 2. (Top) tachyon potential as a function of time, Eq. H5UI for a = 0.1 and 6 = 
(e = 1). (Middle) amplitudes and (bottom) propagation speeds of magnetic field (upper 
three lines) and electric field (lower three lines) for kTo =0.1 (dashed line), kTo = 1 (solid 
line), and kTo = 10 (dotted line). 
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the electric field E = E y (t, x)y is then given by 

E y {t,x) = E(t)e lk[x - WEW \ (53) 

with 



E = \ W B B2 + k~ 2 B 2 , (54) 

W E = W B ~ k- 1 tan" 1 -, . (55) 

\W B B J 

With this ansatz, Eq. I|44|l reduces to 

B-k 2 (w B - e~ 2 V 2 ^j 5 = 0, (56) 

2BW B + BW B = 0. (57) 



Eq. 1(57)1 implies 



^ = i(f)\ (»> 

where Bq is a constant. From Eq. 1)55(1 we also find 

^ = £(§) 2 . (59, 

Inserting Eq. 1)58)1 into Eq. 1)56)1 . we obtain a second-order nonlinear equation 
for amplitude of the magnetic field 

4 

-V 2 



B ^ 



Bo 
B 



B = 0. (60) 



With the knowledge of V(T(f)), like the one given in Eq. (|50|l . we can solve 
these equations for B(t) and -E(i). For V(T(t)) of Eq. JSUJ), the numerical 
solution is shown in Fig. [21 Let us first examine a few limiting cases. 

On the top of the potential V = 1, we get a normal plane- wave solution 
Wb = We = 1/e with amplitudes of B = Bq and £ = £?o/e as expected 
from the original wave equations 1)44)1 1)45)1 . At the top of the potential, the 
equation of states 1)48)1 gives w om = 1/3 and the trace of the energy- momentum 
tensor T en f M up to quadratic terms ()4"§)) vanishes. 

Once the tachyon starts rolling from the top of the potential, the fluctua- 
tions of the gauge fields starts deviating from the normal plane wave solution. 
With the initial condition V(t = 0) = 1 for which 4a 2 = 6, V(T(t)) in Eq. (gOJ 
reduces to 

V(T(t))= 1 (61) 

'1 + fcsinh 2 t/T 
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Assuming that b = (e 2 — l)/e 2 is small, which means the tachyon rolls slowly 
initially, we can solve Eqs. I|60(l and l|58|l to the first order in b: 



B{t) = B Q 



1 



6fc 2 T 2 



W s(t) = - + 77 



t - 



■ sinh 2 t/To 



■sinh 2t/T 



(62) 



2(l + fc 2 T 2 ) 

where we rescaled Bo so that B(0) — Bo- This solution is valid as long as 
6e*/ T ° <C 1. Thus the amplitude initially grows quadratically as the tachyon 
starts rolling. On the other hand, electric field can be calculated from (|54H 
and J5JJ, 

' (i 



E(t) = B 
t 



1 - 



w E {t) 



4(1 
b 



PT 2 ) 



2k i T£ 



8e(l + k 2 Tg) 



[2fc 2 T 2 t 



+ fc 2 T 2 sinh^ t/To) 
fc 2 T 2 )sinh 2t/T ] 



(63) 



Note that E(t) and We tends to decrease in contrast with the behavior of 
magnetic field. This can be confirmed in Fig. [21 As time elapses, w em in 
Eq. H48|) starts decreasing but T on f up to quadratic terms in Eq. (|49|l develops 
nonvanishing contributions. 

At late-time, the value of the potential l|5U|l decays exponentially to zero 

In this case, it is more convenient to work with the original equation (|44H 
from which we see that the time dependence is exponentially suppressed at 
late time. Therefore the magnetic field linearly diverges in general, 

B{t) = B 1 t + B 2 + O{e~ t / T0 ). (65) 

Then the propagation speed decreases according to Eq. 1 (58) 1 (See also Fig.|2J|. 
Since B linearly increases, our quadratic approximation of the Born-Infeld 
action is not valid as t gets large. However, one may still extract some relevant 
physics as the tachyon rolls to the vacuum, V — > 0. The electromagnetic 
waves become frozen and stop propagating. For the magnetic field, short 
wavelengths freeze earlier than long wavelengths. For the electric field in this 
limit, we have from Eq. (|52|) . 

E{t) = Bi/k, (66) 

which is a constant. Note that in (|59|l We contains extra V 2 compared with 
Wb and hence the electric field is frozen more quickly than the magnetic field. 

The density p em (|46(l diverges exponentially from the electric component 
but the pressure p em lE} approaches a finite value, so that w em in Eq. l(4*H) 
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becomes zero. Similarly, the trace of the energy-momentum tensor <|49[) also 
diverges exponentially. Since the Hamiltonian density TLt for the pure tachyon 
is almost a constant, the energy density (|46[) can be compared with that in 
the materials 

^(D-E + H-B). (67) 

Then this material may be approximated by a dielectric material at least 
formally 

D "W' H " B ' (68) 
where its time-dependent dielectric function diverges as time elapses, e{t) w 
l/V(T) 2 ^ oo. 

4 Cosmological Aspects of Rolling Tachyon 

An attractable to pic of the rolling tachyon is its application to cosmol- 
ogy | f f | l^ | f ^ | 14 | 15 | ag i nc Ji ca ted from the beginning The simplest set for 
the cosmology of the rolling tachyon is composed of the graviton and the 
tachyon T. Therefore, the action is read from Eq. 

S = J d^x^—g (J^R - T 3 V(T)y/T+g^%T%Tj . (69) 

Without any computation, we read a few intriguing subjects from charac- 
ters of the tachyon potential in Eqs. ijfJJl 10 (Refer also to Fig. QJ. At initial 
stage of the rolling tachyon, there exists a cosmological constant V(T wO)»l 
so that one expects easily solutions of inflationary universe. At late-time of 
T — > oo, we have tiny nonvanishing but monotonically-decreasing cosmologi- 
cal constant at each instant, which lets us consider it as a possible source of 
quintessence. Bad news for comprising a realistic cosmological model may be 
absence of reheating and difficulty of density perturbation at late-time due 
to disappearance of perturbative degrees reflecting nonexistence of stationary 
vacuum point in the monotonically-decreasing tachyon potential. 

For cosmological solutions, we try a spatially-homogeneous but time- 
dependent solution 

ds 2 = -dt 2 + a 2 {t)dn 2 k , T — T(t), (70) 

where dVi 2 , corresponds, at least locally, to the metric of S* 3 , E 3 , or iJ 3 ac- 
cording to the value of k = 1,0, — 1, respectively. Then the Einstein equations 
are summarized as 

a 2 k k 2 T 3 V(T) , N 

- + - = — 71 
a z a z 3 y/i _ f2 
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a _ k 2 T 3 V(T) / 3 - 2 



T , (72) 



a 3 1 Vj _ y2 \ 2 
and the tachyon equation becomes 

3^ + 1^=0, (73) 

where is equivalent to conservation equation of the energy-momentum tensor. 

For the flat universe of k = 0, cosmological evolution is rather simple 
since both sides of Eq. I|71|l are positive semi-definite. The second equation 
11721 tells us that the expanding universe is accelerating at onset of T(i) « 
since the right-hand side of Eq. (|72|l is also positive definite. As time t goes, T 
exceeds \/2/3 so that expansion of the universe slows down and then finally 

the scale factor a(t) will halt as a(t) constant. This makes the second 
term of the tachyon equation (|73|l vanish at infinite time so that we confirm 
f -> and f -> 1 as T -> oo. 



5 String Cosmology with Rolling Tachyon 

In this section, we study the role of rolling tachyons in the cosmological model 
with dilatonic gravity. In the string frame, flat space solutions of both initial- 
stage and late-time will be obtained in closed form. In the Einstein frame, we 
will show that every expanding solution in flat space is decelerating. 



5. 1 String frame 

We begin with a cosmological model induced from string theory, which is 
confined on a D3-brane and includes graviton g^ Vi dilaton and tachyon T. 
In the string frame, the effective action of the bosonic sector of the D3-brane 
system is given by 

S = At / d^xj~g e~ 2 * (R + 4V„*V#) 
2k 2 J 

-T 3 Jd A x e~* V(T) y f-det{g liV + d ll Td u T) , (74) 

where we turned off an Abelian gauge field on the D3-brane and antisym- 
metric tensor fields of second rank both on the brane and in the bulk. 
Equations of motion for the metric, the dilaton, and the tachyon derived from 
the action (|74^l are 

~ \g„uR = « 2 (T*„ + e*T%), (75) 
4 (V^V^ 1 * - V^^V^ 1 *) + R = k 2 T 3 e* V(T)y/l + V Q TV Q T, (76) 
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l + V a TV«T V ^ v<p VdT ~^ ^ 7 > 



where energy-momentum tensor of the dilaton is 

T* v = ~ [V„V„$ + (V Q $V Q $ - V Q V a $)] , (78) 
and that of the tachyon 

l£ = T 3 ^ =4= = [V M TV,T - + V Q TV a r)] . (79) 



For cosmological solutions in the string frame, we try in addition to the 
metric ansatz l|7(J[l 

$ = (80) 
From the action (|74() . we obtain the following equations 



d 2 fc 



a . 



3 U + ^- 2 l 3 ^" $ J = ' tePT! (81) 

2^ + ^ + ^ - 2 f $ + 2^$ - $ 2 ") = - K 2 e V, (82) 
a a 1 a z \ a J 

4 f$ + 3^>-$ 2 ') -b(° + *+JL) =« 2 eV, (83) 



a 



1 



. 3- - $ T + = =0, (84) 

where tachyon energy density pt and pressure pt defined by T T ^ U = 
diag(-pT,PT,PT,Pr) are 

p T = T 3 -2=!= and p T = "Tg^v/T^T 2 , (85) 
V 1 - T 2 

which formally coincide with Eq. (|2(J|) and Eq. (|22l) . The tachyon equation 
(|84|l is equivalent to the following conservation equation 

Pt + {3H - <S>)T 2 p T = 0, (86) 

where H — a/a is the Hubble parameter. In the absence of detailed knowledge 
of V(T), we will examine characters of the solutions based on the simplicity 
of tachyon equation of state 

Pt = wtPt, wt — T — 1, (87) 
which is exactly the same as Eq. I)2H[I. 
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By defining the shifted dilaton = 2$ — 3 In a, we rewrite the equations 
©-(BU as 

2 -3H 2 + e4 = 2n 2 eiJp T , (88) 



2 

2(H - H0) + 4-^- = K 2 e^ai PT , (89) 
a z 

<j) 2 - 20 + 3H 2 + 6^ = -K 2 ei a ip T , (90) 

Note that y/—g, or a 3 is not a scalar quantity even in flat spatial geometry, 
the shifted dilaton cf> is not a scalar field in 3+1 dimensions. The conservation 
equation l|86|l becomes 

PT + \{'iH -<P)T 2 PT = 0. (92) 

Now Eqs. lj%%|) - lj§?Hl and Eq. lj%7| are summarized by the following two equa- 
tions 

.... fc 

2cj)-(f> 2 + 2H(f)-3H 2 -2H=-lQ—, (93) 

cr 

w T 4> 2 + 4H(j> - 3w T H 2 -4H= -(8 + 6w T )4- (94) 

Let us consider only the flat metric (k = 0) in the rest part of the paper. If 
we express the dilaton ^ as a function of the scale factor a(t), <fi(t) — <fi(a(t)), 
we can introduce a new variable ip such as 

^ = a<j>'=^, (95) 

where the prime denotes the differentiation with respect to a, and the second 
equality shows that ip is the ratio between 4> and H . Then Eqs. (|93|l and (|94ll 
are combined into a single first-order differential equation for ip : 

Aaijj' + {4> 2 -3){w T ^ + 2-w T ) =0. (96) 

From now on we look for the solutions of Eq. (|96|) . Above all one may 
easily find a constant solution ip = =F\/3 which is consistent with Eqs. (|88|) ~ 
(f^TOfl only when px = : 

a(t) =a (lTV3H t)^/^, H(t) = - H ° (97) 

1 =F VoH t 

*(*) = *o - Ml T V3Hot), (98) 
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where H = H(t = 0), ao = a(t = 0), and $ = ®(t = 0) throughout this 
section. However, exactly-vanishing tachyon density pt — from Eq. (|88|l re- 
stricts strictly the validity range of this particular solution to that of vanishing 
tachyon potential, V(T) — 0, which leads to T — oo. The tachyon equation 
(19 111 forces T — and T — 1 so that the tachyon decouples (wt = Pt = 0). 
Therefore, the obtained solution 197|) - (|98|l corresponds to that of string cos- 
mology of the graviton and the dilaton before stabilization but without the 
tachyon. 

Since it is difficult to solve Eq. (|96|) with dynamical wt, let us assume 
that wt is time- independent (or equivalently a(i)-independent). We can think 
of the cases where the constant wt can be a good approximation. From the 
tachyon potential (jSJ, the first case is onset of tachyon rolling around the 
maximum point and the second case is late-time rolling at large T region. In 
fact we can demonstrate that those two cases are the only possibility as far 
as no singularity evolves. 

When wt is a nonzero constant, Eq. (|96|) allows a particular solution 

i> = ^l s fi. (99) 

WT 

This provides a consistent solution of Eqs. H88 |l -19() (l 

-l 



2(2w T - 1) ,„ ^ , w T + 2 



m = *o + ^ T ^ln(l + ^ r H t). (101) 



From Eq. I|88|l and Eq. I|92|) . the tachyon energy density pt is given by 

2(l + mT ) 2 

2 - 2wt - Wt 9 / Wt + 2 \ -t+ 2 

Since the obtained solution is a constant solution of ip, it has only three 
initially- undetermined constants. Specifically, the solution should satisfy $ = 
[(2wt — ^)/wt]H so that the initial conditions also satisfy a relation <J> = 
[(2wt — ^)/wt\Hq. Once we assume general solutions of a(t)-dependent ip with 
keeping constant nonzero wt, they should be classified by four independent 
parameters (00,^0,^0)^0) instead of three in Eq. (|101f> . 

According to the aforementioned condition for valid wt values, the ob- 
tained solution in Eq. (f 1 () If) may be physically relevant as the onset so- 
lution of wt = — 1 (tp = 3). In this case, pT(t) is reduced to a con- 
stant pr(t) = 3e~*°i/Q/K 2 . Comparing this with the definition of pt 
in Eq. I|85() . the initial Hubble parameter Ho is related to the dilaton as 
H = ±ne®°/ 2 y / T 3 /3. Then, with T(t) — 0, the tachyon equation of motion 
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is automatically satisfied and hence Eq. I|1U1|I becomes an exact solution of 
the whole set of equations of motion i|81|) - (|84[l . Since the tachyon field re- 
mains as constant at the maximum of the potential, this solution describes 
the expanding or shrinking solution depending on the initial Hubble param- 
eter, with a constant vacuum energy corresponding to brane tension due to 
tachyon sitting at the unstable equilibrium point. b 

In order to study the behavior of the tachyon rolling down from the top 
of the potential, now we slightly perturb this solution, i.e., look for a solution 
with nonzero but small T dependence. So we treat T as a small expansion 
parameter and work up to the first-order in T. Since the unperturbed solution 
satisfies 3H = $, the tachyon equation of motion 184f> becomes, to the first- 
order in T, 

f-\%. 

This can easily be integrated to 

1 • T 2 

-T 2 = - In V + const = h const, (104) 

2 8 In 2 ' y ' 

where we used the form of the potential near the origin (JSJ . Given the initial 
condition that T = Tq and T = at t = 0, we can solve this equation (jl04[l 
and obtain 

T(t) = T cosh at, (105) 

where a = l/2\An2. Therefore tachyon starts to roll down the potential as a 
hyperbolic cosine function. Taking derivative, we find T = aTosinhat. The 
range for which T remains small is then t < t r = 2%/In2 sinh _ \2v^n2/T ), 
during which the approximation wt ~ — 1 is good. Unless the initial value 
To is fine-tuned, the tachyon follows the onset solution (|101fl - <|102|) f° r t < t r 
and enters into rolling mode. 

For more general solutions, the first-order differential equation H96|l can 
be integrated to 



C 



ip 2 - 3 ftp - V3 N 



75 



(i>-(3) 2 U + \/3, 



13^-3 

(106) 



where C is an integration constant. Note that this algebraic equation does 
not provide a closed form of tp in terms of the scale factor a(t) except for a 
few cases, e.g., w T = 0,-1/(^3- 1/2), -2/(3V§ - 1). 



'The interpretation as expanding or shrinking solution needs to be more careful, since we 
are working in the string frame. Actually the behaviors are reversed in the Einstein frame 
as we will see in the next subsection. 
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Fortunately, for the late-time case of vanishing wt, we can obtain the 
solution in closed form 



V' = \/3 



-V3/2 



c 4 



-a/3/2 ' 



(107) 



Then the scale factor a and the dilaton $ are explicitly expressed as functions 
of time t by solving the equations 189|) - (|90|l : 



a(t) = a 



C+t + 2 



1/V3 



H(t) = 



(C_i + 2)(C+i + 2)' 



$(t) = $ + In 



) (C_£ + 2)(^- 1 )/ 2 
'(C+t + 2)(v / 3+ 1 )/ 2 



(108) 



(109) 



where C± 
Eq. (JHH1) 



(3 =F ^/3)Hq — 2$o. We also read the tachyon density px from 



pr = C+C~e 



(C , _i + 2)(^ / 3+l)/2' 



(110) 



Note that C± should have the same sign from the positivity of the tachyon 
density I|11U|) . Let us first consider that both C+ and C_ are positive. When 
C_ > C+ or equivalently Hq > 0, the scale factor a is growing but saturates 
to a finite value such as a(oo) = ao(C_/C+) 1 / v ^ in the string frame. When 
C_ < C+, it decreases. When C_ = C+, Ho = so that the scale factor is 
a constant, a(t) = ao. For all of the cases, the dilaton $ approaches negative 
infinity. Note that wt — means late-time, the tachyon density decreases 
like pT~l/iast— + oo. Consistency check by using Eq. (|92|l or equivalently 
by Eq. (|91|l provides us the expected result, T — » 1. If both C + and C_ 
are negative, there appears a singularity at finite time irrespective of relative 
magnitude of C + and C_. 



5.2 Einstein frame 

In the previous subsection, it was possible to obtain the cosmological solutions 
analytically in a few simple but physically meaningful limiting cases. To study 
the physical implications of what we found, however, we need to work in the 
Einstein frame. In this subsection, we will convert the cosmological solutions 
obtained in the string frame to those in the Einstein frame and discuss the 
physical behaviors. 
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First of all, let us rewrite the action l|74|) in the Einstein frame of which 
metric g^ u is related to the string metric by a conformal transformation 

Note that we abbreviate the superscript E for convenience in what follows. 
Then the action l(74() is changed to 



S =2 



^2 J d 4 a:\/=<? (i? - 2V M W M $) 
-T 3 f (Px^f^j e 3 * V(T)y/l + e- J *V^rVT. (112) 



Equations of motion are 



1 



R„u ~ \g^R = k\T% + e 3 *Tjj, (113) 
y VM<e = !^ T3 e ^ V{T) 3 + 2e- 2$ V a rV a T 

^ i + e -2*v Q rv«T y <ir" ' 

(115) 

where the corresponding energy-momentum tensors are 

T*„ = ^ (2V„3V„* - ffAII ,V Q $V Q <i>) , (116) 



V{T) 



(117) 

For cosmological solutions, we use again Eq. 1)70(1 and Eq. 1)80(1 and then 
the equations of motion 1114(1 1(115|) become 

" 2 1 K \ P <, + e™ PT ) = \& + ^T z e^V{T) 1 , , (118) 

3 3 y 1 - e -2* T 2 

2 

a 6 

= - J$ 2 + ^T 3 e^VjT)' (119) 
3 6 ^1 - e- 2 *T 2 

$ + 3^4> = -^e 3 *(p T - 2 PT ) = -^T 3 e 3 *X/(T) A=^=2L, (120) 
a 2 2 y 1 - e- 2 *T 2 
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Tachyon energy density pr and pressure px in the Einstein frame arc obtained 
by the replacement f s = e^f in Eq. {gj, 

and p T = ~T 3 V(T)Vl - e'^T 2 , (122) 



Pt/pt = er^f 2 - 1. (123) 

Alternatively, we can obtain the same equations from Eqs. (|81f) - (|84|l by 
using the metric of the form c 

ds 2 = e^i-dt 2 + a 2 (t)dn 2 k ), (124) 

and hence the time t and the scale factor a are related to those in the string 
frame as 

a s = ae*, dt s = e^dt. (125) 

Then the Einstein equations for the flat case (k = 0) in the string frame 
l|5T )l -l(51 )l are converted to 

ff 2 = ^ 2 + ^ 2 e 3 V, (126) 

H = -¥ - \ 2 e *p T T 2 . (127) 

Demanding constant wt is nothing but asking a strong proportionality condi- 
tion between the dilaton and tachyon, T oc e*. Note that the pressure px as 
shown in Eq. (|85(l is always negative irrespective of both specific form of the 
tachyon potential (V(T) > 0) and the value of the kinetic term (e _2 *T 2 < 1), 
and the value of wt interpolates smoothly between —1 and 0. 

First we observe that the right-hand side of Eq. I|126l) is always positive, 
which means that the Hubble parameter H(t) is either positive definite or 
negative definite for all t and it cannot change the sign in the Einstein frame. 
Obviously it is a natural consequence of the weak energy condition. Let us 
first consider the case of positive Hubble parameter, H{t) > 0. Eq. (|127|l 
shows H consists of two terms both of which are negative definite for all t. 
Since H > by assumption, the only consistent behavior of H in this case 
is that H vanishes as t — > oo, which, in turn implies that $ and e* ptT 2 go 
to zero, separately. It also means that H should be a regular function for 

c In this subsection all the quantities are in the Einstein frame except the variables with 
subscript s which denote the quantities in the string frame. 



V{T) 

PT = ^3 , 

VI - e- 2 *T 2 
and thereby wt is given as 

wt 
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all t. In order to find the large t behavior of H(t), one has to study e~*T 
in large t limit which appears in the definition of wt in the Einstein frame. 
Knowing that the functions are regular, it is not difficult to show that the only 
possible behavior is e~*T — > 1 as t — ► oo after some straightforward analysis 
of Eqs. l|12GH|127p . Combining it with the fact that $ and e^prT 2 vanish, 
we can immediately conclude from Eq. 1|126[) that H{t) should go to zero in 
large t limit. 

The asymptotic behavior of fields in case of the positive Hubble parameter 
can be found from the solution l|ll)9|l since wt is essentially zero for large t as 
we just have seen above. The only thing to do is to transform the expressions 
in the string frame to those in the Einstein frame, using the relation Q125|) . 
Therefore, for large i, we find 



a(t s ) = a a (t,)e-*<*'> ~ ^a s0 e-*°(C + t s + 2)^(C_t s 



n/3- 



2) 



/ = / ,lLe~* ~ 2e~*° / dt s {C+U + ■ (128) 

' (C-t s + 2)(^-i)/2 y ' 



One can also identify the initial Hubble parameter Ho in terms of C± as 



H = \e 



1 -^) c - + ( 1 + 4 lr ' 



(129) 



Note that C± have the same sign as the Hubble parameter H . Now, with 
C± > 0, one can easily confirm from the above equation (|128l) that all the 
functions indeed behave regularly. In t s — > oo limit, a ~ t s and t ~ t 2 s so that 
the asymptotic behavior of the scale factor becomes a ~ t 1 ! 2 . This power 
law expansion in flat space is contrasted with the result of Einstein gravity 
without the dilaton where ultimately the scale factor ceases to increase, 
limt—,00 a(t) — > constant. The behavior of tachyon density px can be read from 
Eq. (|110fl with t replaced by t s , which shows that pr ~ t^ 1 ! 2 . Since wt also 
goes to zero, the fluid of condensed tachyon becomes pressureless. Differently 
from ordinary scalar matter where matter domination of pressureless gas is 
achieved for the minimum kinetic energy (T — > 0), it is attained for the 
maximum value of time dependence (e _ *T — * 1 as T — * oo) for the tachyon 
potential given in Eq. JSJl. 

When the Hubble parameter H is negative, the situation is a bit more 
complicated. Since H < always, H becomes more and more negative and 
there is a possibility that eventually H diverges to negative infinity at some 
finite time. Indeed, it turns out that all solutions in this case develop a 
singularity at some finite time at which H — > — oo and a — > 0. These big 
crunch solutions may not describe viable universes in the sense of observed 
cosmological data. Depending on initial conditions, the dilaton $ diverges 
to either oo or — oo and T goes to either oo or zero with the factor e~*T 
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remaining finite. It is rather tedious and not much illuminating to show this 
explicitly, so here we will just content ourselves to present a simple argument 
to understand the behavior. Since the tachyon field T rolls down from the 
maximum of the potential to the minimum at infinite T, it is physically clear 
that T s — e~*T would eventually go to one unless there is a singularity at 
some finite time. Suppose that there appeared no singularity until some long 
time had passed so that e~*T approached to one sufficiently closely. Then 
Eq. i|128|) should be a good approximate solution in this case. However, we 
know that both C± are negative when H < and Eq. 1|128[) is clearly singular 
in this case. We have also verified the singular behavior for various initial 
conditions using numerical analysis. 

As mentioned in the previous section, the tachyon T is decoupled when 
e~*T = l and T = oo. In this decoupling limit, characters of the Einstein 
equations l|126() ~ (|127[l that H' 2 > and H < do not change so that all the 
previous arguments can be applied. Well-known cosmological solution of the 
dilaton gravity before stabilization of the dilaton is 

a(t) = a (l + 3flo*) 1/3 , H{t) ''" 



l + 3.ffo*' 



$(t) 



3£r t), 



(130) 
(131) 



where the (±) sign in Eq. (|131fl is due to the reflection symmetry ($ <-> — $) 
in the equations i|126[) ~ (|127[) . This solution can also be obtained throughout 
a transformation l|125|) from Eqs. (|97l) - (|98|l . For H < 0, it is a big crunch 
solution (a — * 0) which encounters singularity (H — > oo, $ — * ^foo) as t — > 
1/S\Hq\. For H > 0, it is an expanding but decelerating solution. Since 
a ~ i 1 / 3 , the power of expansion rate is increased from 1/3 to 1/2 by the 
tachyonic effect as expected. 

So far we discussed generic properties and asymptotic behaviors of so- 
lutions in the Einstein frame. Now we consider the behavior at the onset. 
The solution (| 1 1|> obtained by assuming constant wt is transformed to the 
Einstein frame as 



a(t) = a 



=>*(«) 



<i>(j 



2(1 -w T ) 

, (wt-2) 2 
2(1 -w T ) 



H t 



2(2u, T - 



H t 



(132) 



where the initial Hubble parameter Hq is related to that in the string frame by 



H 



"HosO- ~ Wt)/wt- Note that Ho and Hq s have opposite signs since 



wt < 0. Therefore the expanding (shrinking) solution in the string frame 
corresponds to the shrinking (expanding) solution in the Einstein frame. For 
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the onset solution with wt = — 1, the tachyon energy density px is a constant 
as before, pr(t) = 3e _3 *° Hq j 4k 2 . Then the initial Hubble parameter is given 
by Hq = ±2Ke 3 *°/ 2 -\/T3/3, which describes the exact solution that tachyon 
remains at the origin as explained in section 2. Under a small perturbation 
tachyon starts rolling down according to Eq. pt)5(l with t replaced by t s . The 
rest of the discussion on the rolling behavior is the same as in the string frame 
and the details will not be repeated here. 

In conclusion the cosmological solution can be classified into two cate- 
gories depending on the value of the Hubble parameter H (t) in the Einstein 
frame. When the initial Hubble parameter Hq is positive, the solution is reg- 
ular and the universe is expanding but decelerating as a(t) ~ \ft while e*W 
vanishes. When Hq is negative, there appears a singularity at some finite time 
t at which the universe shrinks to zero. 

6 Summary 

In this paper we have discussed basic field theoretic properties of rolling 
tachyon and Abelian gauge field, and then cosmological solutions. We have 
shown that both tachyon and Abelian gauge field cannot propagate in terms 
of perturbative modes when the homogeneous tachyon rolled to the infinity 
which corresponds to the minimum of its potential. Throughout the study of 
homogeneous and isotropic universes in string frame, we found initial-stage 
and late-time solutions in closed form in addition to the known solution of 
the gravitons and dilaton in the decoupling limit of the tachyon. We have 
also provided a description of cosmological solutions in the Einstein frame : 
Once the universe starts expanding, then it continues expanding eternally but 
decelerating. 

We conclude with the list of intriguing questions for further study. The 
nonlocal effective action we dealt with is not valid for large value of the tachyon 
field, so genuine string theory description is needed in such regime, including 
tower of massive string spectra. Since both the tachyon and Abelian gauge 
field cannot support propagating physical degrees in a perturbative way, an 
intriguing question may be existence of extended objects in the background 
of the rolling tachyon configuration. In relation with tachyon cosmology, we 
may raise a number of questions. Stabilization mechanism of dilaton should be 
understood and its positive effect to our unsatisfactory decelerating universe 
solutions is carefully investigated. Various topics for the rolling tachyons 
asked in the Einstein gravity should be addressed again in the context of 
dilaton gravity such as existence of inflationary era, possibility as a source 
of quintessence, reheating without oscillating tachyon modes, cosmological 
perturbation and structure formation, and so on. When the Abelian gauge 
field in concerned, one can ask whether or not radiation dominated era is 
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possible. 
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